#!/bin/bash
#PBS -N test_mpi
#PBS -l nodes=8:ppn=20
#PBS -q batch
#PBS -j oe
#PBS -k oe
#cd $PBS_O_WORKDIR

Nuclei=He4
A=4
Z=2

Nshell=11
Lambda=0.0
HbarOmega_Input=20

Int_Type=opt
#/~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~/#
Int_Form=1 #1: Morten's Form 2: Jisp16 Form
data_path="/home/public/heff_v2b_nnlo-opt_bare_cens/"
data_int="nnlo-opt_bare_N11_hw20.dat"
data_orb="nnlo-opt_bare_N11_hw20_orb.dat"
#/~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~/#

#/~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~/#
#Int_Form=2 #1: Morten's Form 2: Jisp16 Form
#data_path="/home/public/heff_v2b_nnlo-opt_cutoff500_coulomb_bare/"
#data_int="v2b_n2lo-opt_coulomb_bare_"${Nshell}"shl_hw"${HbarOmega_Input}
#data_orb=${Nshell}"jorbits.dat"
#/~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~/#

orbfile=${data_path}${data_orb}
intfile=${data_path}${data_int}

out_path="data/result_"
out_data=${out_path}${Nuclei}_N${Nshell}_hw${HbarOmega_Input}.txt

./RPA.exe <<EOF >${out_data}&
${orbfile}
${intfile}
${nuclei}
${Nshell}
${Lambda}
${HbarOmega_Input}
${Int_Form}
${Int_Type}
EOF
